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Abstract 



The energy and electric field dependence of pulse shape discrimination in liquid xenon have been measured in a 10 gm 



two-phase xenon time projection chamber. We have demonstrated the use of the pulse shape and charge-to-light ratio 
simultaneously to obtain a leakage below that achievable by either discriminant alone. A Monte Carlo is used to show 
that the dominant fluctuation in the pulse shape quantity is statistical in nature, and project the performance of these 
^jPj techniques in larger detectors. Although the performance is generally weak at low energies relevant to elastic WIMP 
^ recoil searches, the pulse shape can be used in probing for higher energy inelastic WIMP recoils. 
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1. Introduction 

Two-phase xenon time projection chambers have set 
competitive limits in the search for weakly interacting mas- 
sive particles or WIMPs [3 [5]. These detectors have a 
liquid target and a gas region above with several pho- 
tomultiplier tubes (PMTs) for detecting scintillation and 
electrodes biased at high voltages for establishing charge- 
collecting electric fields. The typical event consists of 
two bursts of light: the primary scintillation ("SI") and 
the proportional (electroluminescence) scintillation ("S2") 
from the drifting electrons exciting the gas [3]. Back- 
grounds can be removed in several ways. With 3D position 
reconstruction capability, events at the edges of the detec- 
tor and multiple scatter events can be removed, eliminat- 
ing most events from external gammas and neutrons and 
surface contaminants. The charge-to-light ratio of back- 
ground electron recoils (from gammas and betas) is higher 
than that of nuclear recoils from WIMPs, allowing for the 
rejection of bulk betas and single Compton-scatters. This 
method has a discrimination efficiency of >99% down to 
2 keVr (recoil energy) [TJ [3]. The measurement of the 
SI pulse shape also provides discrimination between recoil 
types, although generally at a lower efficiency than that 
afforded by the charge-to-light ratio. Dark matter detec- 
tors using pulse shape discrimination (PSD) with Nal [5], 
high pressure xenon gas 0, single-phase liquid xenon [7l|H], 
neon [9], and argon [10], and two phase-argon [TP, 12] have 
been employed or are currently under development, but 



thus far there has been no exploration of PSD in two-phase 
xenon time projection chambers. This paper explores the 
energy and electric field dependence of PSD in such a de- 
tector at energies relevant to dark matter searches. Also 
explored is the use of the pulse shape quantity in conjunc- 
tion with the charge-to-light ratio to achieve a background 
rejection efficiency beyond that achievable by either dis- 
criminant alone. A pulse shape Monte Carlo is used to 
examine the fluctuations of the pulse shape discriminant 
and determine the efficiency of PSD in larger detectors, 
where the distribution in the time-of-flight of photons be- 
comes important. 

2. Scintillation Mechanism 

In condensed noble elements, both excitation and ion- 
ization produced by particle interactions eventually result 
in scintillation by the following processes: 

1) De-excitation 
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The excimer (excited dimer) Rj is presumed to exist in ei- 
ther the singlet or triplet state with measured decay times 
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Figure 1: Diagram of the prototype with two PMTs and five wire 
grid electrodes. The hquid surface sits midway between the gate and 
anode, which are separated by 4 mm. The distance between the gate 
and cathode is 1 cm. 



of 4 ns and 21 ns, respectively [13]. While the decay times 
have been shown to be independent of the stopping power 
dE/dx, the ratio of singlet to triplet states /1//3 has been 
observed to be positively correlated with dE/dx [IHIIII]. 
This provides a mechanism for PSD. The /1//3 of elec- 
tron, alpha, and fission fragment interactions (ordered by 
increasing dE/dx) in liquid xenon have been measured to 
be 0.05, 0.45, and 1.6, respectively [Tl[Til[TS]. For liq- 
uid xenon, a slow tail with a 40 ns fall time is present in 
electron recoils at zero electric field and disappears upon 
application of a field, which strongly suggests that the slow 
tail is due to charge recombination P3] • The tail is absent 
in scintillation pulses of alpha and nuclear recoils, presum- 
ably because the recombination times in dense tracks are 
much lower than the excimer decay times [12 E]. The 
conversion of singlets to triplets (the lower energy state) 
by free electrons has been proposed as the mechanism for 
the positive correlation between dE/dx and /1//3 [13] - 
events with lower dE/dx have electrons that recombine on 
longer time scales and consequently have a greater chance 
of converting states. In this view, the longer recombina- 
tion times associated with lower dE/dx could be shaping 
the pulse in two ways - through the singlet-to-triplet ratio 
and the recombination time. Our measurement probes the 
pulse shape dependence on dE/dx (by observing events of 
different energy) for a particular type of recoil as opposed 
to observing the correlation across different types of par- 
ticle interactions as previous studies have done. 



3. Experimental Setup 



respectively. This results in an SI signal size of ~1 pho- 
toelectron (pe) per keVr for nuclear recoils and ^5 pe per 
keV for electron recoils from 122 keV gammas at zero elec- 
tric field. (The light yield is roughly field independent for 
nuclear recoils and is very field dependent for electron re- 
coils [T7]). The PMT in the liquid collects more primary 
light than the other as most photons undergo internal re- 
flection at the liquid surface. The wire grids are paral- 
lel stretched wires soldered onto copper-patterned Cirlex 
boards. The three lower grids in the liquid are made up 
of 40 micron beryllium copper wires. The top two grids in 
the gas consist of 120 micron gold-plated aluminum wires. 
All grids have a wire spacing of 2 mm. The active volume 
contains 10 g of xenon and has a height of 1 cm and a 
diameter of 3.73 cm. PTFE reflectors (>95% reflectivity 
at ~175 nm [T^) line the walls of the flducial volume to 
boost the light signal and increase its spatial uniformity. 

By careful selection of materials, extensive cleaning of 
parts, baking and pumping of the detector, and recircula- 
tion of xenon through the detector and a hot metal getter 
purifler, we achieved an electron drift length of 30 cm (cor- 
responding to a charge loss of <3%) that remained con- 
stant throughout the experiment. The detector is housed 
in a vacuum-insulated, liquid-nitrogen-cooled cold finger 
cryostat, which maintained a xenon temperature stable to 
within 0.1 K over several months. Three parallel plate 
capacitors built onto the Cirlex rings monitor the liquid 
level, which was stable to 40 microns. 

The signals from the PMTs are passed through a 5x 
SRS SR445A amplifier and fan- in/out, and digitized at 
500 MHz at a resolution of 8 bits. Each PMT signal is 
digitized by two channels (one fine and one coarse) in order 
to extend the dynamic range from single photoelectrons 
to the saturation point of the amplifier. The PMTs were 
measured to give a linear response in SI and S2 over this 
range. The DAQ triggers on SI for events down to ~5 keVr 
and on S2 for lower energy events with a threshold of ~7 
extracted electrons. An energy threshold of 2 keVr is set 
in the analysis by a cut on SI. 

The data was taken with the liquid xenon temperature 
at 187.50 K with the liquid surface sitting midway be- 
tween the anode and gate electrodes (Figure [1}. The field 
above the liquid was kept at 10 kV/cm, which was nec- 
essary for full liquid-to-gas electron extraction [19^. At 
each of five different drift fields - 0.06, 0.52, 0.88, 1.95, 
and 3.96 kV/cm - three types of data sets were taken - 
electron recoil, nuclear recoil, and calibration - using the 
external sources ^^'^Ba, ^^^Cf, and ^^Co, respectively. The 
high energy gammas from the ^^^Cf were shielded with 8" 
of lead. A data set was also taken to calibrate the signal 
in terms of photoelectrons. 



The prototype detector (Figure [ij consists of five wire 
grids and two Hamamatsu PMTs, one of which sits in the 
liquid (R9288) and the other in the gas (R6041). The 
PMTs have quantum efficiencies of 16.8% and 6.0% at 175 
nm, and collect 50% and 10% of the primary scintillation. 



4. Pulse Shape Analysis 

Data processing begins with the identification of scintil- 
lation pulses. Various pulse shape quantities such as the 
pulse height, integral, rise time, and width arc calculated 



2 




■Nuclear Recoil 



0.1 




40 60 
Energy (keVr) 



80 



Electron Recoil 

Nuclear Recoil 



^i|,'2f^>^^^^i.V*. — ^— 



0.1 



20 40 60 

Energy (keVr) 



80 



Figure 2: Scatter plots of fp versus recoil energy at 0.06 kV/cm. 
Top - Nuclear recoils events from the ^^^Cf data set. The electron 
recoil events from this data set have been removed with a cut in 
log]^g(S'2/51). Bottom - electron recoil events from the ^"^^Ba data 
set with the leakage events highlighted. The solid and dashed lines 
represent the electron and nuclear recoil means, respectively. 



and used to classify each pulse as SI, S2 or noise. Cuts are 
applied to remove all but single-scatter events (one SI and 
one S2). The integrals of the SI and S2 pulses are corrected 
for depth dependence and calibrated with 122 keV gammas 
from ^^Co. The field dependence of the charge and light 
yields of these events have been measured previously 
The nuclear recoil energy is Er — Ee/Ceff ■ Se/Sr where 
E^ is the electron recoil energy scaled linearly from the SI 
peak of 122 keV gammas; Se and 5*^ are the electric field 
dependent scintillation yields of electron and nuclear re- 
coils, respectively, relative to their values at zero field jTZ] : 
and Ceff is the nuclear recoil scintillation yield relative 
to that of the 122 keV electron recoils at zero field. The 
energy dependence of Cef / has been measured by several 
groups [ini [211 1211 1231 m] and there is not yet a con- 
sensus on the correct value. For simplicity, we assume a 



constant value of Ceff — 0.2. 

A preliminary study of the PSD performance was done 
with various pulse shape quantities such as the width at 
half-max, integral/height, prompt/total (fraction of light 
at the head of the pulse) , and a quantity called "multi-bin" 
which is a generalization of the prompt/total fraction to 
more than two bins (as described in |25j). The last two 
quantities give equally the best performances. We report 
results using prompt/total as it requires much simpler cal- 
culations. The prompt/total is defined as 



tn-^twindoLi 



S{t)dt 



Qs{t)dt 



(1) 



where S{t) is the signal, tyaindow is the prompt window size, 
ti = to — 50 ns, and tf = to + 300 ns, with to as the time at 
which the pulse reaches 10% of the maximum amplitude 
on the leading edge. The optimal t^jindow depends on the 
electronics and size of the detector, and is 10 ns for the 
prototype. Figure [2] shows a scatter plot of /p as a function 
of recoil energy for the nuclear recoils from the ^^^Cf data 
set. The second scatter plot shows the electron recoils from 
^•^^Ba data set with the leakage events in the nuclear recoil 
acceptance window highlighted. This window is defined 
as the region above the nuclear recoil fp mean, giving a 
nuclear recoil acceptance of 50%. 

The Gaussian-extrapolated leakage is the fraction of 
electron recoil events in the nuclear recoil acceptance win- 
dow as determined by the Gaussian width of the electron 
recoil fp and distance between the electron and nuclear re- 
coil fp means. For all fields, the leakage is calculated for 
discrete energy bins of edges 2, 6, 10, 20, 30, 40, 50, 60, 70, 
80, and 90 keVr. The electron and nuclear recoil energy- 
dependent fp means are characterized by lOth-order poly- 
nomials (these are represented by the lines in Figures [2]) . 
These polynomials are calculated by iteratively fitting to 
the fp means of many overlapping energy bins, subtract- 
ing out the mean, and repeating until the mean-subtracted 
data has no more energy dependence. Fitting to the flat- 
tened bands provides a more accurate measurement of the 
band width since a shift in the mean of fp within an en- 
ergy bin would otherwise lead to an overestimate of the 
band width. Before fitting to the nuclear recoils, a cut 
in charge-to-light ratio logio(S'2/5'l) is applied to remove 
the electron recoil events. The band separation is the sep- 
aration according to the lOth-order polynomials plus the 
remaining offsets as determined by the Gaussian fits to the 
flattened bands. The same leakage calculation procedure 
is used for log;^Q(5'2/S'l). For this discriminant, the nu- 
clear recoil acceptance window is below the \ogio{S2/Sl) 



5. Pulse Shape Monte Carlo 

A pulse shape Monte Carlo was built to gain a better 
understanding of the width of the fp distribution and to 
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Figure 3: Average pulse (left) and probability distribution functions 
(right) for 78 keVr electron and nuclear recoils at 0.06 kV/cm. The 
slight hump following the peak is an electronics glitch of unknown 
origin. 
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Figure 4: Gaussian-extrapolated electron recoil leakage as a function 
of recoil energy for the data set taken at a drift electric field of 0.06 
kV/cm using the pulse shape discriminant fp. 



Figure 5: The energy dependence of the separation between electron 
and nuclear recoil fp fit means (A/p) (top) and electron recoil fit sig- 
mas (bottom) at 0.06 kV/cm. The error bars represent the statistical 
uncertainties. 



project the discrimination in larger detectors. To simu- 
late the SI signal, several functions are needed: the pho- 
ton emission (excimer decay) time probability distribution 
function (pdf), the photon time-of- flight pdf, PMT gain 
pdf, and average single photoelectron response waveform. 
The photon time-of-flight pdf was obtained with a home- 
made photon propagation Monte Carlo. The Rayleigh 
scattering and absorption lengths used in the Monte Carlo 
are 30 cm and 1 m for the liquid tMll^T]) and the corre- 
sponding values for gas are extrapolated by density. The 
beryllium copper and gold-plated grids and PTFE walls 
are assumed to have reflectivities of 0.1, 0.2 and 0.95, re- 
spectively. Although the metal reflectivities are not well 
known as they depend greatly on the surface condition, 
they have negligible effects on the final results. The Monte 
Carlo shows that the time-of-flight is instantaneous for our 
purposes and thus is not included for simulating data for 
the prototype. (However, time-of-flight pdfs are needed 
for the larger detectors discussed later.) The PMT gain 
pdf is taken from the photoelectron calibration data set. 



The photon emission pdfs are obtained by deconvolving 
the photoelectron response from the average SI pulses. 
The average pulses are calculated by lining up the pulses 
according to their tg and calculating the mean voltage for 
each time bin. This was done for each 2 keVr bin from 2 to 
90 keVr. The average single photoelectron response wave- 
form was obtained from the photoelectron calibration data 
set. Figure |3] shows the average pulses and corresponding 
emission pdfs for 78 keVr electron and nuclear recoils from 
the data set taken at a drift fleld of 0.06 kV/cm. Initially, 
the only objective of this experiment was to measure the 
discrimination efficiency of the charge-to-light ratio. Con- 
sequently, little effort was made to match impedances for 
internal electrical lines, resulting in an ever-present glitch 
seen in the average pulses of Figure |3] This glitch gives 
rise to a sharp spike in the photon emission time pdf at 
about 25 ns and has been removed. 
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Figure 6: Electron recoil leakage versus electric field between 60 and 
80 keVr using fp . The error bars represent the statistical uncertain- 
ties. 
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Figure 7: The fit mean value of fp between 60 and 80 keVr versus 
electric field for electron and nuclear recoils. The error bars represent 
the statistical uncertainty. The error bars for the electron recoils are 
too small to be seen. 



6. Results 

Figure |4] shows the electron recoil leakage as a function 
of energy for the data set taken with a drift electric field 
of 0.06 kV/cm. This data set demonstrated the best PSD 
performance among the five taken. Note that an electron 
recoil leakage of 0.5 means that there is no discrimination 
as the nuclear recoil acceptance is 0.5. Not shown are 
the leakage values calculated by counting events. These 
values agree within error with the Gaussian-extrapolated 
leakage values, indicating that any anomalous interactions 
that would have otherwise affected the discrimination are 
either weak or nonexistent. (In comparison, the Gaussian- 
extrapolated and counting leakages for log]^Q(S'2/S'l) gen- 
erally fail to match due to various anomalous events such 
as those that suffer from charge loss.) The error bars of 
the leakage include only the statistical fit uncertainties. 
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Figure 8: Mean fp versus the mean log]^g(S'2/51) for electron recoils 
between 60 and 80 keVr at five different electric fields, with both 
fp and log]^g(S2/S'l) generally increasing with drift field. The er- 
ror bars represent the statistical uncertainties. The error bars for 
logiQ{S2/Sl) are too small to be seen. 



The leakage in the prototype is very high below 20 keVr 
but is <0.05 for events >55 keVr. For comparison, the 
mean electron recoil leakage with log]^g(52/S'l) is ^0.005 
below 20 keVr and ^^0.002 between 50 and 80 keVr for the 
window with 50% nuclear recoil acceptance. The mean fp 
of the electron recoil events in the ^^'^Ba and ^'^^Cf data 
sets agree within error at 70 keVr. Systematic shifts in the 
fp of the electron recoils can be as large as the statistical 
uncertainties. 

Figure |5] shows the separation in fp and the actual and 
Monte Carlo electron recoil band widths as a function of 
energy. These plots reveal that the rise in leakage with 
decreasing energy is due to both a convergence of the two 
bands and an increase in the band width. The merging 
of the electron and nuclear recoil bands can be explained 
by the convergence in dE/dx and drop in recombination 
light for electron recoils. (The electronic stopping pow- 
ers of alpha, electron and nuclear recoils are given in ^.) 
The small differences between the actual and Monte Carlo 
widths indicate that any intrinsic contribution to the fluc- 
tuations in fp must be relativity small and that statistical 
fluctuations dominate. 

Figure [6] displays the fleld dependence of the leakage 
between 60 and 80 keVr. The leakage increases as the 
electric field is increased. The band width changes negligi- 
bly with the electric field since the intrinsic component is 
very small and the statistical contribution does not change 
as the number of photoelectrons per unit recoil energy re- 
mains roughly the same. The trend in the leakage is due 
to the change in the mean fp of the electron recoils as seen 
in Figure |7] The mean fp of the electron recoils possesses 
a strong field dependence below ^1 kV/cm which is simi- 
lar to the behavior of the electron recoil charge yield [T7]. 
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Figure 9: Scatter plots of /p versus log]^Q{S'2/51) (both with respect 
to the nuclear recoil mean) between 70 and 80 keVr at 0.06 kV/cm 
with the solid and dashed straight line representing the bounds of 
the logjQ(S2/51) and 2D windows, respectively. The red diamonds 
and magenta circles highlight the leakage events of the log]^Q(S'2/Sl) 
and 2D windows, respectively. The dashed circles represent the 2a- 
ellipses. 



Figure |8] shows the positive correlation between mean fp 
and log]^Q(S'2/S'l) (a measure of the charge yield) between 
60 and 80 keVr. In contrast, the charge yield and pulse 
shape of nuclear recoils both change little with field [TZ] . 
This strongly suggests that charge recombination plays a 
central role in the pulse shape. Although, a positive cor- 
relation is observed for the mean values across different 
electric fields, the correlation for events at an electric field 
is roughly zero. The correlation is either nonexistent or is 
masked by the statistical fluctuations. 

The electron recoil leakage can be reduced without 
decreasing nuclear recoil acceptance by changing from 
the usual nuclear recoil acceptance window defined by 
logio(S'2/S'l) to one defined in logio(S'2/S'l)-/p space ("2D 
window"). For simplicity, the lower bound of the 50% 
acceptance window is defined as a diagonal straight line 
through the center of the nuclear recoil cluster. The opti- 
mal angle of the boundary and corresponding leakage are 
calculated based on the band separations and electron re- 
coil widths in log]^Q(S'2/S'l) and fp with the assumption 
that the two quantities are uncorrelated. Figure |9] is a 
scatter plot of fp versus log]^Q(S'2/S'l) from the -'^'^^Ba and 

Cf data between 70 and 80 keVr at 0.06 kV/cm with 
the leakage events of the two windows highlighted. Figure 
10 shows the Gaussian-extrapolated leakages with the two 
windows in the prototype (as well as the active volume- 
averaged projected performance for the larger detector dis- 
cussed in the next section) . The reduction in the Gaussian 
leakage is ~6% below 20 keVr and ~84% between 60 and 
90 keVr. 
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Figure 10: Gaussian-extrapolated electron recoil leakage versus recoil 
energy at 0.06 kV/cm with the log j^q (52/51) and 2D windows. 



7. Pulse Shape Discrimination Projections 

The results can be used to estimate the discrimination 
performance of a detector with higher light collection. The 
band separation would be the same. The band width at 
energy i? of a detector with a factor x times the primary 
light collection of the prototype is the sum (quadrature) 
of the Monte Carlo width at a; x plus the intrinsic width 
(difference between the actual and Monte Carlo widths) at 
shows the projected leakage at 0.06 kV/cm 
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E. Figure 

for the prototype with two times the current light col- 
lection efficiency. Although increasing the light collection 
efficiency significantly decreases the band width, the leak- 
age below 20 keVr is still rather high because of the small 
and immutable band separation. 

In addition to simulating the photon propagation for 
the prototype, a next generation two-phase xenon detector 
called the Large Underground Xenon Detector, or LUX, 
and a detector with four times the dimensions of LUX 
(LUXx4) were also simulated. Nominal LUX active vol- 
ume dimensions of 49 cm diameter and 54 cm height were 
used. Both detectors are presumed to have the same light 
collection as the prototype. The time-of- flight distribu- 
tions for fifty (r, z) coordinates were obtained. At each 
electric field, 24,000 SI pulses of each recoil type were 
simulated between and 100 keVr at random points in 
the detectors. For these simulated pulses, t window — 36 
ns is used in the definition of fp. Figure 11 shows the 
detector-averaged electron recoil leakage in LUX, LUXx4, 
and the prototype. The projections for LUX and LUXx4 
use the widths as predicted by the Monte Carlo. The PSD 
performance is very weak for LUXx4 at all energies and 
thus PSD would probably never be used in such a large 
detector except for events near the PMTs whose photons, 
to a large degree, strike the PMTs directly. 

Also projected were the leakages in LUX and LUXx4 



using the 2D window (Figure 10 1. The projections use the 
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Figure 11: Gaussian leakage as a function of recoil energy for the 
prototype, the prototype with two times the primary light collec- 
tion, LUX, and LUXx4 at a drift field of 0.06 kV/cm using the 
discriminant fp. 



fits to log2Q(S'2/S'l) from the real data in conjunction with 
the fits to fp from the Monte Carlo data. The reduction 
in leakage for LUX by changing to the 2D window is ~2% 
below 20 keVr and ~55% between 60 and 90 keVr. The 
improvement for LUXx4 (not shown) is negligible at all 
energies. 

The use of PSD would clearly not be productive in a 
large detector searching for elastic WIMP recoils as the 
efficiency of PSD is weakest at low energies where most of 
the events are expected to appear. However, PSD might 
be useful in inelastic recoil searches. Models for inelastic 
interactions predict a spectrum with few or no events at 
low energies and a peak at ^^40 keVr [5H] UHl 130] ■ 



8. Conclusion 

We have studied the energy and electric field dependence 
of PSD in a two-phase xenon time projection chamber and 
have shown that a reduction in leakage could be obtained 
by using a nuclear recoil acceptance window that utilizes 
both the pulse shape and charge-to- light discriminants. A 
pulse shape Monte Carlo shows that the fluctuations in the 
pulse shape quantity are primarily statistical in nature and 
that PSD in larger detectors is severely weakened by fluc- 
tuations in the photon time-of-flight. We also find that 
improvements in light collection are unlikely to meaning- 
fully improve the discrimination efficiency at lower ener- 
gies due to the small and unchangeable separation between 
the electron and nuclear recoil bands. Although the per- 
formance is weak at the energies relevant to searches for 
elastic WIMP recoils, PSD may be useful in looking for 
higher energy inelastic WIMP recoils. 
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